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Abstract 

The rebellious voter model, introduced by Sturm and Swart (2008), is a variation of the 
standard, one-dimensional voter model, in which types that are locally in the minority 
have an advantage. It is related, both through duality and through the evolution of its 
interfaces, to a system of branching annihilating random walks that is believed to belong 
to the 'parity-conservation' universality class. This paper presents numerical data for the 
rebellious voter model and for a closely related one-sided version of the model. Both 
models appear to exhibit a phase transition between noncoexistence and coexistence as 
the advantage for minority types is increased. For the one-sided model (but not for the 
original, two-sided rebellious voter model), it appears that the critical point is exactly a 
half and two important functions of the process are given by simple, explicit formulas, a 
fact for which we have no explanation. 
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1 Introduction 



The study of parity preserving branching-annihilating particle systems has a long history in 
the physics literature. Their main interest lies in the fact that they provide a rare example 
of systems that exhibit a phase transition of survival/extinction type (the physics literature 
usually speaks of a transition between an 'active' and an 'absorbing' phase) that does not 
belong to the directed percolation (DP) universality class. 

Indeed, the phase transition of these models appears to belong to a different 'parity conser- 
vation' (PC) universality class, characterized by its own critical exponents. The first example 
of a phase transition belonging to this universality class was described in [GKT84j ; this con- 
cerned a one-dimensional model for interfaces, motivated by solid state physics. Numerical 
simulations for parity preserving branching-annihilating particle systems on Z d in dimensions 
d = 1,2,3 and also on Sierpinski gaskets can be found in [TT92J. Further simulations for a 
one-dimensional model were made in [ Jen94j . This paper also contains the conjecture that the 
order parameter critical exponent (3 equals 13/14; in another paper the author has conjectured 
that (3 = 1. The latter value was also found in [IT98] , while more recent literature seems to 
agree on a value of (3 ~ 0.92 or slightly higher, see e.g. [HinOO, OS05]. Our present work is 
relevant to the question of determining the value of (3 since our proposed explicit formula for 
the survival probability p of the one-sided rebellious voter model, if it is correct, implies that 
[3 = 1. 

Field-theoretic (renormalization group) work on the PC universality class was carried 
out in [CT98], who found two critical dimensions. The (usual) upper critical dimension is 
d c = 2, above which one should find mean-field exponents. There is, however, a second critical 
dimension d' c ~ 4/3 such that only in dimensions below d' c there is a nontrivial absorbing 
phase. From the field-theoretical perspective, the existence of a second critical dimension 
leads to considerable difficulties, to which the authors of [CCDDM05J reacted by considering 
non-perturbative techniques. 

As usual, mathematical theory lags far behind theoretical physics and occupies itself with 
more basic and partially different questions. In particular, renormalization group techniques 
are generally nonrigorous and the existence of critical exponents or their universality below 
the upper critical dimension are from the mathematical perspective unproved conjectures. 
Nevertheless, in recent years, some rigorous mathematical work has been carried out on mod- 
els (supposedly) belonging to the PC universality class. Sudbury [Sud90] studied a 'double 
branching annihilation random walk' (DBARW), which is a one-dimensional model for parity 
preserving branching-annihilating particles. The DBARW that can be exactly solved, in a 
certain sense, but unfortunately the model is in the absorbing phase for all nontrivial values 
of its parameter, hence there is no 'true' phase transition. (Nevertheless, its solution has been 
used in papers such as [CT98| to draw conclusions that presumably hold for the whole PC 
universality class.) No solvable models seem to be known so far that exhibit a true phase tran- 
sition. (Below, we will present a possible candidate for such a model, the one-sided rebellious 
voter model.) 

Strong impetus to the mathematical study of models from the PC universality class 
came from the paper of Neuhauser and Pacala {NP99] , which proposes a variation on the 
(i-dimensional range-A voter model motivated by applications in population biology. Their 
model is governed by two parameters ao, ot\. Here aj < 1 models the interspecific competition 
rate, i.e., the death rate of organisms of type i due to competition with organisms of the 
other type, while the intraspecific competition rate is supposed to be one. In the symmetric 
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case ao = ct\ = a, their model is dual to a system of parity preserving branching- annihilating 
random walks, with branching rate 1 — a. For d = 1, N = 1, up to a rescaling of time, this dual 
model is the DBARW; results for the latter imply that the model of Neuhauser and Pacala 
with d = 1, N = 1 exhibits noncoexistence for all a > 0. For N > 2 and d = 1, their model is 
supposed to exhibit a nontrivial phase transition between coexistence and noncoexistence. For 
d = 2, their model is supposed to exhibit noncoexistence only for a = 1 while in dimensions 
d > 3 noncoexistence holds for all a, in line with the predictions made in [CT98J. 

Much of this lacks rigorous proof. First of all, it is not rigorously known that coexistence 
is a decreasing property in a, or, in terms of the dual model, that increasing the branching 
rate leads to more particles in the system. Obvious as this may seem, finding a rigorous 
proof seems hard. Second, there exists no rigorous proof of noncoexistence (or equivalently, 
the existence of an absorbing phase for the dual model) except in the trivial cases a = 1 
(pure voter) or d = 1, N = 1 (related to DBARW). On the other hand, there exist a number 
of rigorous results proving coexistence (or equivalently, the existence of an active phase for 
the dual model). Neuhauser and Pacala proved that their model exhibits coexistence for a 
sufficiently close to zero for all values of d and N except d = 1, N = 1. In [CP 06J. it has been 
proved that the model of Neuhauser and Pacala exhibits coexistence in dimensions d > 3 for 
a sufficiently close to one (the case of intermediate a is still open due to the lack of a proof 
of monotonicity in a). The analogue result in d = 2 is proved in [CMP09]. 

Some similar models were introduced in [BEM06J, who studied parity preserving branching- 
annihilating random walks which allow for more particles at one site, and |SS08a| . who in- 
troduced the rebellious voter model that will be studied in the present paper. A simple 
numerical simulation of that model, run for illustrative purposes, seemed to suggest that the 
critical point of this model is exactly a half, which motivated the present work. In fact, our 
results show that the critical value for this particular model is not ^, but a closely related, 
one-sided version of the model seems to have its critical point exactly at a = ^ and may in 
some sense be exactly solvable. 

2 Set-up and basic properties 

2.1 Definition of the models 

Let 

{0, 1} Z := {x = (x(i)) i& : x{i) G {0, 1} Vi € Z} (2.1) 

be the space whose elements x are infinite arrays of zeroes and ones, situated on the integer 
lattice Z. We will be interested in Markov processes taking values in {0, 1} Z . We will mainly 
focus on the rebellious voter model introduced in |iS_S08aj, and a variation on this model, which 
we call the one-sided rebellious voter model. For expository reasons, we will also shortly discuss 
two other models, the disagreement and swapping voter models. 

In the one-sided rebellious voter model, initially, each site i£Z has a type x(i) € {0, 1} 
assigned to it, and these types are updated according to the following stochastic dynamics. 
Independently for each site i E Z, one constructs times < T\{i) < T2(i) < ■■■ according to a 
Poisson process with rate 1, i.e., T\(i) — 0, T2(i) — Ti(i), ^(i) — T2(i), . . . are i.i.d. exponentially 
distributed with mean 1. At each time t = Tfc(i), with k = 1, 2, . . ., the type of site i is updated 
according to the following rules. With probability a, the site i looks at the site i — 1 on its 
left, and if the type of site z — 1 is different from the type of i, then site i changes its type. 
With the remaining probability 1 — a, the site i looks at the two sites i — 2 and i — 1 on its 
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left, and if the type of site z — 2 is different from the type of i — 1, then site % changes its type. 
If we let X t {i) denote the type of site i at time t, then (X t ) t >o is a continuous-time Markov 
process taking values in {0, 1} Z , which we call the one-sided rebellious voter model. Another 
way of describing its dynamics is to say that in this model, for any i £ Z, the process X makes 
the transition 

x i — y with rate 

(2.2) 

«l{i(i-l)/x(!)} + (1 - a)l{a;(i-2)^a:(i-l)}3 

where 1a denotes the indicator function of an event A and for any x G {0, 1} Z and A C Z, we 
let 

x\j): { l ,.; rU ' ! (2.3) 
v ' I x(j) if j # A y 1 

denote the configuration obtained from x by changing the types of all sites in A. 

The original rebellious voter model as introduced in [SS08a] is similar to the process de- 
scribed above, except that sites do not only look to the left, but also to the right when deciding 
whether to update their type. More precisely, in this case, at each time Tfc(i) as defined above, 
the site i decides, with probability 1/2 each, to look either to the left or to the right. If the 
site looks to the left, then its type is updated as before. If the site i looks to the right, then 
with probability a, it looks at the site % + 1 on its right, and if the types of % and i + 1 are 
different, then i changes its type. With the remaining probability 1 — a, the site i looks at 
the two sites i + 1 and i + 2 on its right, and if the types of i + 1 and i + 2 different from each 
other, then i changes its type. Another way of describing these dynamics is to say that in this 
model, for any i £ Z, the process makes the transition 

x \— > x^ with rate 

2 a {l{x(i-l)^x(i)} + l {x(i)^x(i+l)}) ( 2 -4) 

+ 5(1 - OL) (l{a;(i_2)/a;(i-l)} + l{x(i+l)^x(i+2)}) ■ 

We note that our description of the rebellious voter model differs a factor 1/2 in the speed of 
time in comparison to the original definition in [SS08a . 

A third and fourth model, that we also shortly discuss in this introduction, are the dis- 
agreement and swapping voter models. The disagreement voter model evolves as 

x 1— > x^ with rate 

a(l{z(j-l)^r(i)} + + (1 - &)l{x(i-l)^x(i+l)}, 

while the swapping voter model evolves as 

x ^ w ith rate a(l{x(i-i)^z(i)} + l{»(i)/x(i+l)}) > 

x ^, x {i,t+i} with rate (1 - a)l{x(i)^(i+i)}- 

In case of the rebellious voter model and its one-sided counterpart, the parameter < a < 1 
can be interpreted as the (interspecific) competition parameter. For a = 1, the model is a 
standard, one-dimensional voter model, as first introduced in [CS73. IHL75] , For a < 1, one 
can check that the dynamics give an advantage to types that are locally in the minority. As 
a consequence, at least on a heuristic level, decreasing a should make it harder for any type 
to die out. The rebellious voter model has been introduced in |SS08a| in an attempt to model 
the distribution of two closely related species where competition between organisms belonging 
to different species is less strong than competetion between organisms of the same species. 



(2.5) 



(2.6) 
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2.2 Interface and dual models 

The models introduced so far are in two ways related to parity preserving branching annihi- 
lating particle systems. First, if X is any of the models defined in (|2.2[) and (|2.4p ~ (j2.6p . then 
setting 

Yt(i) ■= l { x t (i)?x t (i+i)} (t > 0, % € Z) (2.7) 

defines a Markov process taking values in {0, 1} Z that we call the interface model associated 
with X. (The sites i such that Xt(i) ^ Xt(i+ 1) are called interfaces, or also kinks or domain 
walls.) Indeed, in case of the one-sided rebellious voter model, it is not hard to see that Y 
jumps as 

y ^ y {i,i+i} W ith ra te al {y(i)=1} + {1 - a)l {y(i _ 1)=1} , (2.8) 
where we use notation defined in fj2.3|) . For the two-sided process, the dynamics of Y are given 

y i — y with rate 

( 2 - 9 ) 

5a(l{»(i)=l} + l{y(i+l)=l}) + s(l - «)(l{y(i-l)=l} + l{y(i+2)=l}) 

(see [SS08al Section 1.2]), while for the disagreement voter model we get the 'swapping anni- 
hilating random walk' (SARW) 

y^y{hi+i} with rate a(l {y(i)=1} + l{j,(»+i)=i}) + (1 - a ) l {y(i)^y{i+i)} > ( 2 - 10 ) 
and for swapping voter model we get 



y ^ yW+i} with rate a (l {y{i)=1} + l{ J/ (*+i)=i}] 

y y{i-hi+l} with rate (1 - a)l {y{i)=1} . 



(2.11) 



Following [Sud90| . we call this latter model the double branching annihilating random walk 
(DBARW); very similar models (usually in discrete time) have been called BARW2 or BAW 
in the physics literature (see, e.g., |TT92| ). 

In each of these interface models, since Y always flips the types of two sites at once, it is 
easy to see that Y preserves parity, i.e., if Y is started in an initial state Yq which contains a 
finite, even (resp. odd) number of ones, then Yf contains an even (resp. odd) number of ones 
for each t > 0. In particular, if Y is started with an odd number of ones, then the ones can 
never completely die out. 

There is a second relation between our models and parity preserving branching-annihilating 
particle systems, namely, through duality. Recall that two Markov processes X and Y with 
state spaces Sx and Sy are dual to each other if there exists a function ip defined on Sx x Sy , 
called duality function, such that 

W[^{X u y)\ = E»[^(x,y t )] (t > 0), (2.12) 

where (resp. M y ) denotes expectation with respect to the law of the process X (resp. 1") 
started in Xq = x (resp. Yq = y). A necessary and sufficient condition for (|2.12p is that 

Gx4>(-,y)(x)=GYi>(x,-)(y) (xeS x ,yeS Y ), (2.13) 

where Gx (resp. Gy) denotes the generator of A (resp. Y). In typical applications, one consid- 
ers tp such that the functions - ,y)) y€Sy and (tp(x, ■ )) x€Sx are distribution determining. 
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In this case, the transition probabilities of X and Y determine each other uniquely through 

In our case, each of the voter models X we have introduced is a cancellative spin system 
in the sense of |Gri79| . hence it is know that for each of these models, there exists a Markov 
process Y which is also a cancellative spin system and which is dual to X with the duality 
function ip(x,y) := (— l)l 22/ l, where for any x,y,z € {0, 1} Z we let xy(i) := x(i)y(i) denote the 
componentwise product of x and y, and we write \z\ := ^2iz(i). In this case, (I2.12|) implies 
that 

P[|X t Y | is odd] = F[\X Y t \ is odd] (t > 0) (2.14) 

whenever X and Y are independent (with arbitrary initial laws), where in order for our 
expressions to be well-defined we must assume that either \Xq\ or |1q| is finite. 

A special property of the rebellious voter model (that motivated its introduction in [SS08a]) 
is that its interface model and dual model coincide, i.e., if X is the rebellious voter model, 
then (I2.14p holds for the same Y, with dynamics described in (12, 9h . that also describes the 
interfaces of X. One way of checking this is to verify (|2,13p ; an alternative proof, based on the 
graphical representation of cancellative spin systems, is described at |SS08al formula (1.8)]. 
For the one-sided rebellious voter model, the dual and interface models do not coincide, but 
its dual is the mirror image of its interface model, i.e., the one-sided rebellious voter model 
satisfies (|2.14p for a model Y whose dynamics are described by (compare (|2.8p ) 

y^yb-^t} with rate al {y ^ =1} + (1 - a)l {y ( i+V)=1} . (2.15) 

It turns out that the dual (in the sense of cancellative spin systems) of the disagreement model 
is the DBARW (which is the interface model of the swapping voter model) while the dual of 
the swapping voter model is the SARW (which is the interface model of the disagreement 
voter model) (see |SS08at Section 2.1]). Note that in the SARW there is no branching, which 
has far-reaching consequences for both the disagreement and swapping voter models. 

2.3 Coexistence, survival and interface tightness 

It is easy to see that the constant configurations and 1 are traps for the various voter models 
we have introduced. By definition, we say that a Markov process taking values in {0, 1} Z 
exhibits coexistence if there exists an invariant law that is concentrated on configurations that 
are not identically zero or one. It is known that nearest-neighbor voter models on Z d exhibit 
coexistence if and only if d > 2; in particular, for a = 1, our models, being one-dimensional 
voter models, do not exhibit coexistence. On the other hand, our models exhibit coexistence 
for a = since in this case one can check that product measure with intensity a half is an 
invariant law. 

The question is whether there is a nontrivial phase transition between coexistence and 
noncoexistence. For the disagreement and swapping voter models it turns out that this is not 
the case. Indeed, it can be rig orously proved |Sud90l INP991 ISS08aj that these model exhibits 
noncoexistence for all a > 0. On the other hand, it appears that the rebellious voter model 
and its one-sided counterpart do exhibit a nontrivial phase transition between coexistence and 
noncoexistence at some < a c < 1. 

We will mainly be interested in two functions of our processes. It is said that a voter 
model survives if there is a positive probability that the process started with a single one 
never gets trapped in 0. In view of this, for the (one- or two-sided) rebellious voter model 
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with competition parameter a, we will be interested in the survival probability 

p{a) := F 5a [X t ^ Vt > 0] , (2.16) 

where P 5 ° denotes the law of the process started in Xq = So, with 5o(i) := l{j=o}- It seems 
plausible that for many models of the type we are considering, survival is equivalent to co- 
existence. Indeed, using the fact that its interface and dual models coincide, this can be 
rigorously proved for the rebellious voter model [SS08a, Lemma 2]; the same proof works for 
the one-sided process. One of the main aims of the present paper is to obtain good numerical 
data for the function p. 

The function p quantifies, in a sense, how strongly the process exhibits coexistence (see 
also formulas (12. 7j) and (13, ip below). We need a similar formula to quantify noncoexistence. 
For any i £ Z, let £ {0, 1} Z , defined as 

denote a 'Heaviside' configuration of zeros and ones that are 'completely separated'. Following 
[CD95], we say that a model X (with given parameter a) exhibits interface tightness if 

X (a) := lim F x ° \X t = xf for some % € Zl > 0, (2.18) 

<->oo 

i.e., if the process started in a Heaviside configuration spends a positive fraction of its time in 
such configurations. It is intuitively plausible that interface tightness implies noncoexistence. 
It is rigorously known that the swapping voter model exhibits interface tightness for all a > 
[SSOSbj . 

The main aim of the present paper is to obtain good numerical data for the functions p 
and Xi both for the rebellious voter model and for its one-sided analogue. 



3 Main results 
3.1 Methods 

In our simulations, we start the interface process Y of either the rebellious voter model or 
its one-sided counterpart with an odd number of ones on an interval of N sites with periodic 
boundary conditions. Note that in this case, because the number of ones is odd and because of 
the periodic boundary conditions, we can no longer represent Y as in (12. 7j) . but the dynamics 
in (|2.8p and (|2.9p . as well as the duality relation (|2. 14[) . still make sense. Since we start with 
one particle, because of parity conservation, the system cannot die out. Letting Yx, denote the 
process in equilibrium, we assume that for large N, the following approximations are valid: 

(i) p(a) « 2P[y oo (0) = 1] =|E[|r 0O |], 

(ii) x(a)^F[\Y 00 \ = l]. 

The arguments why formula (|3.ip (ii) should hold are obvious from (|2.7p and (|2.18p . The 
justification of formula (|3.ip (i) is more subtle. Let Y 1 / 2 denote the interface process of the 
rebellious voter model on Z, started in an initial law such that the (Yo(*))iez are independent 
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Figure 1: The interface model Y of the rebellious voter model for a = 0.4, a = 0.5, a = 0.51 
and a = 0.6, started with a single interface. Space is plotted horizontally and time vertically, 
with a black spot at a point (i, t) indicating that Yt(i) = 1. Total time elapsed in each picture 
is 1800. The horizontal bar is 500 sites long. 



with P[Yo(*) = 0] = P[YoW = 1] = I' anc ^ ^ De rebellious voter model started in 
X = 5 Q . Then by duality (l2TT4j) . 

P[y// 2 (0) = 1] = P[|y t 1/2 X | is odd] = P[|y 1/2 X t | is odd] = \¥[X t ± 0] , (3.2) 

1/2 

where denotes the constant zero configuration. It is known that the law of Y t converges 
as t — > oo to an invariant law, which is in fact the unique spatially homogeneous, nontrivial 
invariant of the process [SS08a, formula (1.12) and Theorem 3 (b)]. Taking the limit t — > oo 
in (|3,2p . we see that for this invariant law, formula (|3.ip (i) holds with equality. Since it seems 
reasonable that the invariant laws of our finite systems approximate the invariant law of the 
infinite system, this justifies (|3. 1 1) (i). Similar arguments apply to the one-sided model. 

3.2 A first qualitative look 

Before we present our numerical data for the functions p and x from (|3.ip . we first take a look 
at the qualitative behavior of the interface model Y of the rebellious voter model and its one- 
sided analogue. In Figure [IJ we have graphically plotted Y for the two-sided process, started 
with a single one, with time running upwards and a black spot at a point (i, t) indicating that 
Yt(i) = 1. The figures shows the process for four values of a, where a c ~ 0.51 is the estimated 
critial value for this process (see Section I3~3l below) . 

The pictures show that for a < a c , all ones are contained in an interval with edges that 
grow at an approximately linear speed, and that inside this interval, ones occur at some 
approximately constant equilibrium density. On the other hand, for a > a c , the process 
spends a positive fraction of its time in states where there is just a single site in state one, 
which indicates that the rebellious voter model exhibits interface tightness (see (I2.18P ). What 
happens exactly at a = a c is less clear from these pictures, but our numerical data, presented 
in Sections 13.31 and l4~3l below, suggest that at a = a c , the edge speed is zero and the rebellious 
voter model exhibits neither coexistence nor interface tightness. 
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Figure 2: The interface model of the one-sided rebellious voter model for a = 0.3, a = 0.5 
and a = 0.6. Total time elapsed in each picture is 1800. The horizontal bar is 500 sites long. 



The picture for the one-sided process is very similar, except that now, due to the nature of 
the process, the process as a whole tends to drift to the right at some approximately constant 
speed; see Figure EJ 

3.3 The two-sided model 

We simulated the interface model Y of the rebellious voter model on an interval of N sites 
with periodic boundary conditions by slowly increasing or decreasing a from some initial value 
«b to some final value a e during a time interval of length T. We then divided our time interval 
into n equal pieces and plotted the average value of 2|Y|/iV (for the function p) or the fraction 
of the time that \Y\ = 1 (for the function x) against the average value of a, for each of the n 
time intervals. 

This method has both advantages and disadvantages. One of its main advantages is that 
it allows one to quickly obtain data for a very large number of values of the parameter, which 
is useful if one is interested in curve fitting or in estimating derivative functions such as in 
Figure [TBI The method introduces obvious errors due to the nonzero speed of varying the 
parameter, but it is not difficult to get a rough idea of the size of these effects. In fact, one 
can get a good idea of the size of the errors just by looking at the curves (see Section 13.51 for 
a more detailed discussion) . A disadvantage of the method is that the size of the errors varies 
a lot along the curves since relaxation times increase as one approaches the critical point. 
One can partially compensate for this by doing detailed simulations near the critical point. 
(A better approach, which we have not pursued, would be to vary a at a suitable, precisely 
chosen nonconstant speed.) 

The statistical errors in our data for the function p appear to be approximately normally 
distributed (with a variance that depends strongly on a), but those in the data for x have a 
skew distribution that is very far from a normal one (this is somewhat visible from Figure [7]) . 
This non-Gaussian behavior is not unexpected. Indeed, if we start the interface model Y in a 
state with |1q| = 3 particles, then it is known that for the model with a = 1, 



P[inf{s > : \Y S \ = 1} > t) ~ t~ 3/2 



as t — > oo. 



(3.3) 
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Figure 3: The functions p and x for the rebellious voter model. Data obtained with parameters 
as in (15H). 
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Figure 4: The function p(a) — (1 — c±a)/(l — C2«) for the values c\ = 1.957, c 2 = 0.973 (left), 
ci = 1.958, c 2 = 0.975 (middle) and ci = 1.959, c 2 = 0.977 (right). 



It is seems natural to conjecture that this is also true for a < 1, hence the process started 
with a single particle makes excursions into states with more than one particle that have a 
duration which has a finite first moment but infinite second moment. 

In Figure [3] we have plotted approximations for the functions p and x obtained by our 
methods, using the parameters: 



for p : N = 2 15 , n = 2 9 , T = 10 8 , a h = 0, a e = 0.55, 
for x ■ N = 2 15 , n = 2 9 , T = 10 12 , a b = 0.99, a e = 0.49. 



(3.4) 



The values p(0) = 1 and x(l) = 1 are expected on theoretical grounds. Indeed, it is easy 
to check that the rebellious voter model with a = never dies out (recall (|2.16p ) while in the 
pure voter case a = 1, the process started in a Heaviside state stays in such a state for all 
time (recall (12381) ). 

The numerical data for p are well fitted by a linear fractional function of the form 

1 — C\OL 



p(a) 



with ci = 1.958 ± 0.001 and c 2 = 0.975 ± 0.002, 



(3.5) 



1 — C2a 

(see Figure H]). Assuming this sort of function fitting is justified, we arrive at an estimate for 
the critical point of 

a c = 0.5107 ±0.0003. (3.6) 
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Figure 5: Estimation of the critical point of the rebellious voter model. Above: The function 
p(a)/((Xo — a) for the values ao = 0.508 (left), ao = 0.510 (middle) and ao = 0.512 (right). 
Below: the function x(a)/(a — ao) for the values ao = 0.508 (left), ao = 0.510 (middle) and 
ao = 0.512 (right). Our simulations for p use the parameters N = 2 15 , n = 2 7 , T = 10 9 , 
ab = 0.55, a e = 0.45 and those for \ use the parameters in (|3.4p . 



This estimate is to be viewed with caution, since this is based on extrapolation of data for 
(approximately) < a < 0.35, assuming that the linear fractional form (|3.5|) holds for all 
a < a c . A more robust method (see Figure [5]) based on our best data for the functions p and 
X near the critical point leads to the estimate: 

a c = 0.510 ±0.002. (3.7) 

(See also Figure [9j) It seems that linear fractional functions do not fit the numerical data for 
X really well. 

Our data strongly suggest that p(a c ) = = x( a c), which implies that at criticality, the 
process exhibits neither coexistence nor interface stability. Summarizing, the picture that 
emerges from our numerical data is as follows: 

There exists a critical value a c ~ 0.510 such that the rebellious voter model exhibits 
coexistence if and only if a < a c and interface tightness if and only if a > a c . The 
function p from (I2.16D has a finite negative slope and is concave on [0,a c ], and is 
approximately given by a linear fractional function of the form (j3.5[) . The function 
X from (I2.18D has a finite positive slope and is concave on [a c , 0]. 

3.4 The one-sided model 

We have run the same simulation as described in the previous section, with parameters as 
in (|3.4p . also for the one-sided rebellious voter model. The resulting approximations for the 
functions p and x are plotted in Figure [6l 

At first sight, Figure [6] looks extremely similar to Figure [31 A closer inspection reveals, 
however, that the two plots are not identical. Indeed, it seems that for the one-sided model, 
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Figure 6: The functions p and \ f° r the one-sided rebellious voter model. Data obtained with 
parameters as in (|3.4p . 



the functions p and x are described by the explicit formulas: 

p( a ) = v 1 - 2a and = V (2 ). (3.8) 

1 — a a 

In particular, one has the symmetry p(l — a) = x( a ) an d the critical parameter seems to be 
given by 

q c = 0.500 ± 0.002 = i (?) (3.9) 

In Figure we have plotted the differences between the explicit functions in f|3.8|) and our 
data for p and Xi respectively. The systematic deviations of p and x from the proposed 
formulas near the critical point probably stem from finite size effects (compare Figures [11] and 
[T2l below). A more detailed comparison of p with its proposed explicit formula is shown in 
Figure [H Because of the small difference between the left and right edge speed near the critical 
point (see Section [4. 3D . these detailed simulations are sensitive to the direction in which a is 
varied (compare Figure [T4l and the discussion in Section l3.5p . In particular, the systematic 
deviation between our proposed explicit formula and the curve produced by lowering a in 
the left plot of Figure [8] seems to stem from such effects. This systematic deviation almost 
disappears after giving the system more time to relax (see the right plot of Figure E|) . 

We guessed our formula for p after estimating the first few derivatives at a = and the 
formula for x by analogy with p. We do not know of any theoretical reason to expect these 
formulas. A detailed comparison of the two-sided and one-sided model near the critical point 
is shown in Figure [9] 

Our proposed explicit formula for p, if it is correct, implies that the order parameter 
critical exponent (3 for our model is one. Since this is different from the recently published 
values = 0.92(2) (HinOOj and j3 = 0.95(1) [OM06j for the PC universality class, we have 
used our data to obtain a direct estimate for the value of f3. A plot of log(/o) as a function of 
log(a c — a) yields an approximate linear graph as a approximates a c from below, the slope of 
which should be /3. To estimate this slope, in Figure [10] we have plotted log(p) — (3 log(o c — a) 
as a function of — log(a c — a), both for the two-sided and one-sided model, for different values 
of j3. These plots are suggestive of a value of (3 that is perhaps closer to 0.92 than to 1, but 
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Figure 7: Difference between the explicit formulas in (|3.8p and the functions p and x f° r the 
one-sided rebellious voter model, with detail of the left and right side of the picture. 




0.46 0.48 0.5 0.52 0.54 0.46 0.48 0.5 0.52 0.54 



Figure 8: Plots of p(a) — max{0, (1 — 2a)/(l — a)} for the one-sided rebellious voter model 
near the critical point. All plots use the parameters N = 2 15 , n = 2 7 . On the left: T = 10 9 
and «b = 0.445, a e = 0.55 resp. «b = 0.55 and a c = 0.45. On the right: T = 4 • 10 9 and 
«b = 0.445, a e = 0.55 resp. = 0.55 and a e = 0.445. 
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Figure 9: A detailed simulation near the critical point of the functions p and x f°r the two- 
sided process and the one-sided process. On the left: our simulations for p use the parameters 
N = 2 15 , n = 2 7 , T = 4 • 10 9 , a b = 0.445, a e = 0.55. Our plots for X show detail of a 
simulation with the parameters N = 2 15 , n = 2 9 , T = 10 12 , Qb = 0.99, a c = 0.5. 



this is probably due to the fact that our simulations are not precise enough to get sufficiently 
close to the critical point. In fact, our data become unreliable at a c — a < 0.004, which 
corresponds to the point — log (0.004) 5.5 in Figure [TU1 Note that in Figure [U for the one- 
sided model, the same data were shown to be consistent with our proposed explicit formula 
for p in (|3.8j) . which, if true, implies that /3 = 1. We conclude that our data are not good 
enough to convincingly distinguish between /3 = 0.92 and (3 = 1. 

3.5 Finite size effects 

To gain more insight into how close our numerical data are to the real functions p and x, we 
tested the effect of varying the system size N. In Figure [TT1 we have plotted our approximation 
(see (I3.ip ) of the function p for different values of the system size N. The pictures for the 
two-sided and one-sided process are extremely similar, except for a shift of the critical point. 
In view of Figure [TTJ it seems that we can rule out the possibility that the observed differences 
between the two-sided and one-sided process are due to finite space effects (and a hypothetical 
slower convergence for the two-sided model). 

Note that in Figure [TT1 our approximations to p near the critical point become more 
rough as the system size is increased. This can be understood due to two effects. On the 
one hand, near criticality, the random variable {Y^l assumes the values 1,3,5,7,... with 
approximately equal probabilities (see Section [4.2|) . which means that the number of ones in 
the system fluctuates from being close to one to a positive fraction of all sites in the lattice. 
As the system size gets large, this means huge fluctuations, which are moreover slow since at 
criticality, the edge speed is zero. 

In Figure [T2"| we have plotted our approximations for the function x f°r various values 
of N. To counteract the 'roughening' effect we have just described, in these simulations, we 
have increased the total time T together with the system size N. Again, the pictures for the 
two-sided and one-sided process are similar except for a small shift in the critical point. 



p(oc) (one-sided) 
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Figure 10: Estimation of the order parameter critical exponent (3. Plots of log(p) — (3 log(a c — a) 
as a function of — log(a c — a). For the right choice of f3, the asymptotic slope of these 
functions should be zero as the parameter on the horizontal axis tends to infinity. Above: 
plots for the two-sided model with j3 = 0.92 (left) and f3 = 1 (right), based on the estimate 
a c = 0.51. Below: the same for the one-sided model, based on the estimate a c = 0.5. Each 
figure combines two plots, one based on simulations with the parameters N = 2 
T = 4 • 10 9 , 
alpha («b = 
model). 
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«b = 0.445, a e = 0.55, the other instead using N = 2 and a smaller range of 
0.5, a e = 0.52 for the two-sided model and «b = 0.49, a e = 0.51 for the one-sided 




Figure 11: Effect of the system size on our approximation of the function p from (I3.ip . Plotted 
are our approximations for p using the parameters n = 2 7 , T = 10 s , = 0.55, a c = 0.45, and 
the system sizes N = 2 7 , 2 s , 2 9 , 2 10 , 2 , 2 12 , respectively. On the left: the two-sided model. 
On the right: the one-sided model. 
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Figure 12: Effect of the system size on our approximation of the function \ from f)3. 1 1) . 
Plotted are our approximations for \ using the parameters n = 2 7 , ab = 0.6, a c = 0.49-0.5, 
and the system sizes N = 2 10 , 2 11 , 2 12 , 2 13 , 2 14 and times T = 10 9 , 2 • 10 9 , 4 • 10 9 , 8 • 10 9 , 16 • 10 9 , 
respectively. On the left: the two-sided model. On the right: the one-sided model. 
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Figure 13: Effect of the duration of our simulations on our approximation of the function p 
from (|3,ip for the one-sided model. Plotted are our approximations for p using the parameters 
N = 2 10 , n = 2 7 , o b = 0.6, a e = 0.4, and the total elapsed times T = 10 5 , 10 6 , 10 7 , 10 s , 
respectively. 



We have also run simulations for different values of the total elapsed time, to investigate 
the effect of this parameter on the quality of our data. Since short times mean the system 
does not have enough time for time averages to reach their equilibrium values, the main effect 
of increasing T is to smoothen our approximated functions. For the function p, this effect is 
demonstrated in Figure [T3j A second effect of choosing a too short time T is that the numerical 
functions 'lag behind' in that they show values belonging to an a that lies somewhat in the 
past. This effect is demonstrated in Figure [T4l In this case, our approximations near the 
critical point depend on the direction in which we vary a: the graph produced by increasing a 
'overshoots' the critical point while the graph produced by lowering a picks up too late. Also 
note the 'hook' at the beginning of the graph started with a below the critical point, which 
is due to the fact that we start with a single one and the finite edge speed needs several time 
steps to fill out all space. This sort of effects only occur when T is small relative to N; in all 
of our simulations (except those in Figure I14D we have done our best to choose our parameters 
such that these effects in minimal. In case of doubt, we have run the same simulation in both 
directions for comparison (see Figure [8]). In some plots, we have removed a small 'hook' at 
the beginning. 
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Figure 14: Effects of choosing the time too small compared to the system size, in simulations of 
p for the one-sided model. The two curves were produced by raising resp. lowering a between 
the values ctb resp. a e = 0.4 and a e resp. «b = 0.6. The other simulation parameters were 
n = 2 7 , N = 2 13 , and T = 10 6 . 

4 Other functions of the process 
4.1 Harmonic functions 

We have obtained numerical data for a number of other functions of our processes. We have 
mainly concentrated on the one-sided model, in view of the apparent explicit formulas for the 
functions p and % i n this case. We have not been able, however, to find explicit formulas 
for any other functions than these two. Nevertheless, our additional functions show some 
interesting phenomena, which we discuss here. In the present section, we discuss a harmonic 
function for the (one-sided) rebellious voter model. In the next two sections, we look at the 
equilibrium probabilities of seeing 3,5,7,... ones in the interface model, and edge speeds, 
respectively. 

Recall that if X is a Markov process with state space Sx, then a function / : Sx — >• M is 
called a harmonic function for X if, for any initial state, the process (/(-Xt))t>o is a martingale. 
At least for finite state spaces, this is equivalent to the statement that Gxf = 0, where Gx 
denotes the generator of X. In general, a duality between Markov processes translates the 
invariant laws of one process into harmonic functions of the other process. More precisely, if 
X and Y are dual Markov process, with duality function ip, and p is an invariant measure for 
Y, then using (I2.14|) it is easy to see that setting 



defines a harmonic function / for X. Note that multiplying / with a constant will not change 
the fact that it is harmonic. 

To apply this to the one-sided rebellious voter model, we first need to_introduce some 
notation. In this section, we let Xt denote the rebellious voter model, Xt the one-sided 
rebellious voter model, and Xt the mirror image of the latter, i.e., the model that jumps as 





(4.1) 



x i — y a;™ with rate 




(4.2) 



17 



We let Yi and Yt denote the interface models of Xt and Xf, respectively, i.e., Yt is the process 
with dynamics described in (j2.8[) and Yt is the process in (j2.15|) . Then Yt is dual to Xt and Yt 
is dual to X t , in the sense of (12.141) . 

We consider Y t on an interval_of N sites with periodic boundary conditions, started with 
an odd number of ones, and let Y^ denote the process in equilibrium. For each x € {0, 1}^, 
we define 

? P[|xyoo| is odd] 

fx,N{a) := — -. (4.3 

P[^oo(0) = 1] 

— ► 

By our previous remarks, f x N\°i), as a function of x for fixed a, is a harmonic function for 
the one-sided rebellious voter model Xt, on an^interval of N sites with periodic boundary 
conditions. We will be interested in the limit of f x ,N{ot) as N — > oo, when x remains finite and 
fixed. In view of this, we only write down^only the part of x that is nonzero. For example, we 
write /iioi,jv(oO t° denote the function f x ,N{ot) where x is any element of {0, 1}^ such that 
(x{i),x(i + l),x(i + 2),x(i + 3)) = (1, 1,0, 1) for some i E {0, . . . , N - 1} and x(j) = for all 
j 7^ i, i + 1, i + 2, i + 3. Note that since the law of Y^ is translation invariant (modulo N), 
this definition does not depend on the value of i. 

Our simulations suggest that the functions /a j\r(a) converge as N — > oo to a nontrivial 

limit function f x (a), which, as a function of x for fixed a, is a harmonic function for the 

— ► 

one-sided rebellious voter model Xt on Z, started with finitely many ones. If a < a c , then on 
theoretical grounds one may expect that 

f x {a) = §M where ^(a) := [X t ^ Vt > 0] (4.4) 

— > 

denotes the probability that the one-sided rebellious voter model X t started in the initial state 

x survives. For a > a c , in the regime where interface tightness holds, we can relate f x {a) to 

the invariant law of the process Yt, started with a single one and viewed from its left-most one 

(see formula (|4. 19[) below). Since we expect f x (ot) to be harmonic we expect that Gf x (a) = 0, 
— » 

where G is the generator of the one-sided rebellious voter model, which is given by 

Gf x = ^2 (<xl{x(i-l)jtx(i)} + (1 _ a)l{x(i-2)^(i-l)}) {f x {i) - fx)- (4.5) 
i 

— » — * — > — > 

In particular, the facts that Gf\{a) = and Gfn(a) = lead to the relations 

« (/0(«) " /i(a)) + (/li(«) " /i(«)) + (1 - «)(/ioi(a) - /i(a)) =0, 
.— > -> . . — ► — ► . . — ► -» . 

a (/in (a) - /11(a)) + (/i(a) - /11(a)) + (1 - a)(f lwl (a) - /11(a)) =0, 



(4.6) 



where /0(a) = and /i(a) = 1. An inspection of our numerical data shows that these 
equations are satisfied within the precision of our simulations. 

Our numerical data for the functions f x (a), for some simple choices of the pattern x, are 

— > 

shown in Figure [151 Note that fi(x) (not shown in the figure) is by definition identically one, 

— > 

because of the normalization chosen in (|4.3j) . A surprising fact is that the functions f x (a) 
and their derivatives (shown in Figure [T6j) appear to continue smoothly across the critical 

point. Note that the normalization in (|4.3h is crucial here, since above a c , both the numinator 

— > 

and the denuminator tend to zero as N — > 00. For a = 1 one has f x (a) = \x\, which is a 
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Figure 15: The functions f n , f ni , / mi , . . . (left) and f n , f wl , f 10Ol , . . . (right). The functions 
in the left figure increase to the limit l/p(a) while the functions in the right figure increase to 
2 — p(a). Both limits (theoretical curves) are plotted together with the data. The functions 
/nil fwi are based on data obtained with the parameters N = 2 , n = 2 8 , T = 10 8 , 

2 s , T = 10 9 , a b 
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0, a e = 0.5 (left half of the pictures) and N = 2 
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0.99, 



a e = 0.51 (right half of the pictures). The other functions are a combination of simulations 



using the parameters N = 2 
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Figure 16: Derivatives of the functions from Figure 1151 On the left ^fn(ct), ^fni(o:), ■ ■ ■■ 

On the right ^/n(a), ^7101(0)) ^/iooi(«)) • • •• Figure based on the same data as Figure [T5l 
The derivatives are estimated with a quadratic Savitzky-Golay filter using 11 data points. 
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well-known harmonic function for the pure voter model. For a = 0, one has f x (a) = 1 for all 
nonzero x, reflecting the fact that the process with a = never dies out. 
A few regularities may be observed from Figures [15] and [TH Let 

x n ■= 111111111 and x' n := 10000001. (4.7) 

n ones n-2 zeros 

One expects that for a < a c , 

lim p Xn (a) = 1 and lim p x > (a) = 1 - (1 - pi(a)) 2 . (4.8) 

n— ¥oo n— too 

Hence, in view of (j4.4|) . one expects that 

lim fx n (a) = Pii®)' 1 and lim f x > (a) = 2 - pi(a), (4.9) 

n— ¥oo n— >oo " 

which is indeed what we observe. 

From Figure [TBI we moreover observe that 

i/^(«)L=l = 2(n-l) (n>2). (4.10) 

— > — > — > 

This can be explained on theoretical grounds. Indeed, if we write G = aG vot + (1 — a)G re b, 

— » — > 

where G vo t and G re b are the (one-sided) voter and rebellious parts of the generator, then a 
simple calculation shows that 

G mt ^l(a)\ a=1 = G rch f x (l). (4.11) 
Setting g x := ^/x(a)| a=1 , this yields 

511 + 9% ~ 2pi = G vot gi = G rc b/i(l) = 2, 

ffin + ffi - 2fifn = G vot g n = G rcb /n(l) = 0, 

— > — > — ► 
3iiii + <7n — 2(7m = G vo t3iii = G rc b/in(l) = 0, 

etcetera 



(4.12) 



Since f%{a) = and /i(a) = 1 we have g$ = g± = and we can solve g x inductively for 
x = 11, 111, 1111, . . . from pi5]> . 

The fact that the critical point and various other functions of the invariant law are different 
for the one-sided and two-sided model implies that the law of Y^ cannot be mirror symmetric. 
Indeed, we observe that the functions /ion(o) and /1101(a) are not identical, as show in 
Figure OH 

4.2 Frequencies for three and more particles 

As in the previous section, letting denote the equilibrium interface process for the one-sided 
rebellious voter model, we have plotted in Figure [18] (picture on top left) the functions 

Xk (a) =P[|y oo | = k] (fc = 1,3,5,...)- (4.13) 
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0.2 0.4 0.6 0.8 1 

— ► — > 

Figure 17: The functions /ion and /noi for the one-sided model, demonstrating the asym- 
metry of Yoq. Combination of data obtained with the parameters N = 2 , n = 2 8 , T = 10 8 , 
a b = 0, a e = 0.5 and N = 2 15 , n = 2 8 , T = 10 9 , a b = 0.99, a e = 0.51. 




Figure 18: The frequency = lP[|^oo| = k] of seeing k particles for k = 1, 3, 5, . . . in the 

one-sided model. The second and third plot show the functions 9k(a) := Xk+2{oi) / Xk{o) and 
^fc(a) := #fc+2(a0/#fc(a) for the first few values of k = 1, 3, 5, — The fourth plot shows detail 
of the third plot. Data obtained with the parameters for \ listed in (I3.4p 



21 



In particular, xi is the function x from f)3. 1 1) (ii). Except for the case k = 1 (see (|3.8p ). we have 
not found any simple explicit formulas that fit these curves. Nevertheless, some regularities 
may be noted. In particular, if we define 

6 k (a) := 1—- and <p k (a) := - . , [k = 1, 3, 5, . . .), (4.14) 

Xk(a) 6 h [a) 

then it seems that the functions converge very fast to a limiting function $oo(c0 ■ — 
linifc^oo 9k(a). (See Figure [T8l picture on top right; the functions 63,65,... all seem to fall 
on top of each other.) Likewise, the functions 4>k converge very fast to the function that is 
identically one. (See Figure \TE[ two bottom pictures.) It is hard to obtain sufficiently precise 
data, but our best data (shown here) suggest that the function fa is not identically one, al- 
though it is close. This means that the functions #3 and 65 are probably not identical, even 
though they almost fit each other in the picture on top right. Our numerical data are well 
fitted by the linear /constant functions: 

0i(a)wco + ci(a-i) with c = 0.9115 ± 0.0015 and a = -0.13 ± 0.01, 

2 (4.15) 
03 (a) « 1.0027 ±0.0005. 

and 4>k( a ) = 1 f° r k > 5. We note that if one assumes that the explicit formula for Xi( a ) 
in (13. 8p is correct, and one knows the functions 4>k{ot) for each k = 1, 3, . . ., then using the 
fact that Y^=oX2n+i{ot) = 1 one has enough equations to solve the functions Xk( a ) f° r each 
k = 1,3,.... 

The limiting function 600(a) := limfc_ s . 0O 9k (a) seems to be be strictly decreasing and 
concave on [^,1] and satisfy #oo(|) = 1> #oo(l) = 0. This implies that the distribution of 
the number of particles in equilibrium has an exponentially decaying tail for each a > |. In 
particular, the mean number of particles 

H(a) :=E[|yoo|] (4.16) 

is finite for all a > | and diverges as a I | . 

The picture for the (two-sided) rebellious voter model is extremely similar, except that the 
critical point is no longer 0.5 and we have no explicit formula for Xi( a )- 

4.3 Edge speeds 

Let Yt be the interface model of the (two- or one-sided) rebellious voter model started with a 
finite number of ones and let 

It := inf{i G Z : Y t (i) = 1} and r t := sup{« G Z : Y t (i) = 1} (4.17) 

denote the position of the left-most and right-most ones, respectively. Our simulations suggest 
that there exist constants v-(a) < v + (a), called the left and right edge speed, respectively, 
such that 

It ft 
lim — = v-(a) and Urn — = v+(a) a.s. (4-18) 

t— >oo t t— >oo t 

Set 

Z t (i):=Y t (l t + i) (i>0, t>0). (4.19) 

Then Z% describes the process Yt as seen from the left-most one. It seems reasonable to assume 
that the law of Z% converges as t — > 00 to some equilibrium law. Let Z^ denote the process 
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Figure 19: Edge speeds for the rebellious voter model (left) and its one-sided counterpart 
(right). The pictures show the speed of the right edge and left edge as a function of a. 
Parameters used are N = 2 11 , n = 2 9 , T = 10 9 , Ob = 0, a c = 0.5 (left half of each picture) 
and N = 2 11 , n = 2 9 , T = 10 11 , a h = 0.999, a e = 0.5 (right half of each picture). 



in equilibrium. The left edge speed v-(a) can be obtained as the equilibrium expectation 
E[/(Z OC3 )] of a suitable function / which sums the sizes of all possible changes of It multiplied 
with the rate at which they occur. 

We have simulated the process Z% on a finite interval of ./V sites by neglecting all particles 
that fall off the right edge of the interval due to the dynamics of Yj or the shifts in If. This 
process does not preserve parity, but if N is large enough the probability that all particles 
annihilate each other is very small. For most of our data points, such events never occurred, 
and when they occurred (especially near the critical point) they were so rare that they likely 
had no big influence on our estimate of the speeds. By slowly lowering or raising a in our usual 
fashion and keeping track of all changes in Zt that correspond to a change in It we obtained 
numerical data for the left edge speed v-(a) and, in a similar fashion, also for the right edge 
speed v+(a), both for the two-sided and one-sided model. 

The results are shown in Figure [T9l One has v-(a) < v+(a) if and only if a < a c and the 
functions u_(a) and v+(a) are strictly increasing resp. decreasing on [0, a c ). For the two-sided 
model, by symmetry, v-(a) = — v+(a) and hence V-(a) = = v+(a) for a > a c , while for the 
one-sided model one has < v-(a), but otherwise the two pictures are remarkably similar. 
A detailed simulation near the critical point (shown in Figure [201) shows that the estimates 
for the critical points which one obtains from these simulations are consistent with our earlier 
estimates. Apart from the obvious values v-(0) = and u_(l) = = 1 (for the one-sided 

model) and V- = v + = on [a c , 1] (for the two-sided model) we are not able to say anything 
explicit about the curves v—(a) and v+(a). 

5 Conclusions 

We have simulated two one-dimensional models of parity preserving branching and annihilating 
random walks, which arise as the interface models of the (two-sided) rebellious voter model 
introduced in |SS08aj and the one-sided rebellious voter model introduced in the present 
paper, respectively. Both models appear to go through a phase transition between an active 
and an inactive phase as the parameter a (which is one minus the branching rate) is increased 
from zero to one. In most aspects, the processes behave very similarly both qualitatively and 
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Figure 20: Detail of Figure OS Parameters used are N = 2 15 , n = 2 9 , T = 10 11 , with 
«b = 0.52, a e = 0.5 for the two-sided model and = 0.51, a e = 0.49 for the one-sided model. 
The curves have been smoothed, corresponding to (effectively) n = 2 6 . 



quantitatively. 

A peculiar property of the one-sided model, however, is that both the particle density p 
and the equilibrium probability of finding a single particle Xt defined in ()2. 16j) and (|2.18p . 
appear to be given by explicit formulas, and the critical value appears to be exactly a c = \. 
We have no idea why the one-sided modeLwhich has less symmetry than the two-sided model, 
should be more tractable than the latter HI In fact, most hypothetical explanations that come 
to one's mind should apply to the two-sided model as well, for which our numerical data 
convincingly show that the formulas in (|3,8p do not fit. We have tried without success to find 
explicit formulas for other functions of our processes. 

An important question is whether our proposed formulas in (|3.8p for the functions p and 
X of the one-sided model are exact or only good approximations. In formulas (|3.5p and (|4.15p . 
we have seen examples of explicit formulas containing some 'strange' constants that fit the 
numerical data for certain functions of the one-sided and two-sided model, respectively, within 
the precision of our simulations. In these examples, we believe the given formulas are probably 
not exact but more or less coincidentally close to the real functions. On the other hand, after 
trying for a considerable while to find explicit formulas for other functions of the processes, 
as we did, without success, one really starts to appreciate how well the simple formulas in 
(|3.8p fit our data. The explicit formula for p has moreover been upheld in repeated more 
precise simulations near the critical point. Thus, our present position is that these formulas 
are probably exact, even though we have no theoretical explanation for this. 

An important consequence of our proposed formula for the particle density p of the one- 
sided rebellious voter model, if it is exact, is that it implies that the order parameter critical 
exponent /3 for this model is one. By the principle of universlity, one then expects this 
to be true for the whole PC universality class. It is presently generally believed that (3 is 
somewhat smaller than one. Recent estimates are (3 = 0.92(2) [HinOO] and f3 = 0.95(1) 
[OM06], although older numerical work in [IT98] yielded an estimate consistent with (3 = 1. 

1 Perhaps the only aspect in which the one-sided model is potentially simpler than the two-sided model is 
that in the former, information is passed only from left to right. A third model, which we have only studied 
briefly, in which voter model updates look to the right and rebellious updates look to the left, seems to have a 
critical point a c = 0.510 ± 0.003 and functions p, \ that are close, but not identical to those of the two-sided 
rebellious voter model. 
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A thorough discussion of the reliability of these estimates lies out of the scope of the present 
paper. We stress that while our proposed formula for p strongly suggests that (3 = 1, the most 
convincing data for the correctness of this formula are obtained some distance away from 
the critical point. Our best data near the critical point are not good enough to convincingly 
distinguish between (3 = 0.92 and (3 = 1; see Figure [TUl 

It is currently believed that for the PC universality class (in dimension one), (3/v± = i. If 
also (5 = 1, this suggests that the static critical exponents of this universality class are simple. 
(The dynamical critical exponents, such as those related to the edge speeds, could still be 
nontrivial.) In this context, it is interesting to note the effect, observed in Section [4.11 that 
the harmonic function f x {o>) defined there appears to continue smoothly across the critical 
point, which also suggests simple critical behavior. 

We conclude this paper with some suggestions for further work. First of all, the proposed 
formulas (j3.8|) ask for a theoretical explanation. This does not seem to be easy. The task 
would be easier if one could find explicit formulas for more functions. In particular, if one 
could find explicit formulas for the harmonic functions f x (ct) define in Section 14. 1\ then it 
would be straightforward to prove that these are indeed harmonic, which probably could 
be used to prove (|3.8p as well. For those interested in such problems, we have made the 
data underlying our figures available in the supplementary material. Even in the absence of 
an explicit formula, these harmonic functions (and their observed smooth behavior) seem a 
fruitful object for further theoretical considerations. 

Our paper also calls for a critical evaluation of the evidence collected so far about the order 
parameter critical exponent (3 of the PC universality class. Finally, it raises the question if 
there are more models in this universality class that are (or appear to be) explicitly solvable. 
Our present work suggests that such models might in particular be expected in the class of 
one-sided models, where information is passed in one direction only. 
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